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The role of the interface potential on the effective mass of charge carriers is elucidated in this work. We develop 
a new theoretical formalism using a spatially dependent effective mass that is related to the magnitude 
of the interface potential. Using this formalism we studied Ge quantum dots (QDs) formed by plasma 
enhanced chemical vapour deposition (PECVD) and co-sputtering (sputter). These samples allowed us to 
isolate important consequences arising from differences in the interface potential. We found that for a higher 
interface potential, as in the case of PECVD QDs, there is a larger reduction in the effective mass, which 
increases the confinement energy with respect to the sputter sample. We further understood the action of O 
interface states by comparing our results with Ge QDs grown by molecular beam epitaxy. It is found that the 
O states can suppress the influence of the interface potential. From our theoretical formalism we determine 
the length scale over which the interface potential influences the effective mass. 

PACS numbers: 73.21.La, 78.67.Hc, 68.65.-k, 73.90.+f 


I. INTRODUCTION 

Si and Ge nanostructures (NSs) remain amongst the 
most important materials for photovoltaic and electronic 
applications 1,2 . Si is extensively used in the microelec¬ 
tronics industry partly due to the importance and stabil¬ 
ity of its oxide 3 . The Ge oxide states can create compli¬ 
cations in NSs. Fabrication of Ge NSs in a SiC >2 matrix 
often results in a high concentration of interface defect 
states 4-6 , because Ge can easily occupy the Si position 
within the Si 02 matrix due to its low formation energy. 
Sub-oxide interface states give rise to a defect photolumi¬ 
nescence (PL) band that can dominate quantum confine¬ 
ment (QC) effects 5 . However, with careful control over 
the fabrication' and characterization 8 conditions, Ge NSs 
do clearly demonstrate QC effects. Furthermore, because 
of a larger Bohr radius in Ge compared to Si, Ge has the 
advantage of greater tunability in the gap energy (Eg) 
due to QC 9 . Moreover, it was recently demonstrated 
that interface defect states may enhance the efficiency of 
photovoltaic devices due to preferential hole trapping 9 . 

The nature of its interfacial structure is arguably the 
defining parameter of a NS. Fundamentally, the intro¬ 
duction of an interface modifies bulk parameters 5 and 
introduces new phenomena. Interface states can break 
valley degeneracy 10 . Dangling bonds, stress, and de¬ 
fect states elicit the degree to which the interface po¬ 
tential confines charge carriers 5 . Hydrogen versus oxy¬ 
gen termination at the interface influences the oscilla¬ 
tor strength 11 and creates polarization effects at the 
interface 12 . Nanostructure sensing devices rely on the 
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charge configuration at the interface, which ‘communi¬ 
cates’ with the environment 13 and determines the con¬ 
ductivity within the NS 14 . Nonetheless, the role of the in¬ 
terface remains largely misunderstood 15 . The exact band 
energy alignment at the interface remains a topic of in¬ 
tense research 16 , as it is difficult to establish the exact 
chemical environment at the interface 1 '. Defect engineer¬ 
ing at the interface can have a complicated effect on the 
QC of the charge carriers 18,19 , which can be difficult to 
model theoretically 20,21 . 

The crystal potential changes across the interface of 
a NS and thus the carrier effective mass (EM) changes, 
which is modelled by invoking the Bastard-type bound¬ 
ary conditions 22-25 . There is experimental 6,26-28 and 
theoretical 29,30 evidence that the EM should depend on 
the NS dimension. Understanding how the EM varies in 
low-dimension remains a challenge as experimental mea¬ 
surements are model dependent 31,32 and it is difficult to 
theoretically scale the EM to low dimensions 6,33 . Den¬ 
sity functional theory (DFT) was used to calculate the 
EM from the density of states 29 . However, it is not clear 
how to model the interface nor the excited states within 
DFT 25,34 . 

A common pitfall in theoretical modelling is the lack 
of information regarding the interface atomic potentials, 
which leads to extensive use of H-terminated surfaces at 
the expense of the experimental details 5 . Sub-oxide inter¬ 
faces are typically modelled using OH termination 35,36 . 
The local density approximation is used to correct for 
the local charge environment to obtain agreement with 
the experimental optical gap 5,35 . However, there are re¬ 
cent theoretical results using embedding matrices, such 
as SiC 37 . Ref. 33 embedded Si NSs in a Si02 matrix 
by inserting O atoms between Si-Si bonds followed by a 
randomizing procedure leading to defect and stress free 
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structures. Therefore, relaxed assumptions are required 
to obtain agreement with experiment and the defect land¬ 
scape cannot be readily included 38,39 . 

Owing to the rich variety of phenomena that the inter¬ 
face can give rise to, understanding the contribution of 
the individual components can be difficult. Indeed, the¬ 
oretical studies have suggested 5,6,40 Eq ~ D~ x , where 
1 < x < 2 and D is the NS diameter, which has been 
demonstrated experimentally 5,41,42 . In this work, we elu¬ 
cidate the precise role of the interface potential on the 
EM. An analytic formalism for the EM and interface 30 
was developed and applied to study Ge quantum dots 
(QDs) grown by molecular beam epitaxy (MBE) 43 . The 
parameters developed during our study on MBE Ge QDs 
are transferred to the present study on rf-magnetron co¬ 
sputtering (sputtering), and plasma enhanced chemical 
vapour deposition (PECVD) Ge QDs. The PECVD and 
sputter samples were chosen for this work because they 
were fabricated under similar experimental conditions, 
but exhibit a measurable variation in the chemical struc¬ 
ture of the interface. Both sputtered and PECVD QDs 
have the same dimensional dependence as the MBE QDs, 
which allows us to isolate the role of the interface poten¬ 
tial on the EM. 

Comparing PECVD and sputter samples, it is found 
that as the interface potential increases the EM is re¬ 
duced thus increasing the confinement energy. We deter¬ 
mine the length scale at which the EM is modified due to 
the interface potential. Furthermore, by comparing with 
MBE Ge QDs we find that O interface states can sup¬ 
press this mechanism. These results are important for 
solar cell applications, as we recently demonstrated that 
defects at the interface can be engineered to enhance the 
light harvesting capabilities in Ge QDs 9 . 


II. EXPERIMENT 

In this work, we focus our analysis on experimental 
results from two fabrication methods for Ge quantum 
dots (QDs) that produce different interface potentials. 
Thin films of Ge-rich Si0 2 were deposited using either 
PECVD 8 or sputtering 9 . Thermal annealing was car¬ 
ried out at 800 °C or 600 °C in N 2 for PECVD or sput¬ 
tered samples, respectively, to promote Ge QD nucle- 
ation. In Refs. 8 and 9, the experimental details of diam¬ 
eter control and determination are discussed along with 
crystallinity measurements. Both PECVD and sputtered 
QDs were determined to be in an amorphous state with 
a bulk gap energy {Eq{ oo)) of 0.8 eV. The average diam¬ 
eter in the sputter samples varies from 2 to 4 nm, while 
in PECVD samples it varies from 3.5 to 8.4 nm. Optical 
absorption measurements were performed to assess the 
variation in Eg with QD diameter. Ge NSs embedded 
in an oxide matrix typically demonstrate a PL spectrum 
dominated by a high concentration of sub-oxide related 
defects at the interface between the Ge NS and the oxide 
matrix, thus masking the QC related PL band 5 . In Refs. 


8 and 9, the absorption spectrum was considered over 
the PL spectrum, because defects states in the Ge0 2 in¬ 
terface made it difficult to measure the variation of Eq 
with QD diameter from the emission spectrum 44 , which 
is more sensitive to mid-gap defects 45 . Finally, it is worth 
noting that due to NS size distribution and carrier leak¬ 
age there is some inherent error in associating the NS 
experimental Eq{D) with the correct experimental D , 
see Ref. 5 for more details. 

From the absorption spectra, the oscillator strength, 
and Eg(D) was found in Refs. 8 and 9 for each D using 
the Tauc method. The applicability of the Tauc method 
to NSs is discussed in detail in Ref. 9. The Tauc method 
is based on a semi-empirical model for amorphous sys¬ 
tems, which assumes parabolic band edges and optical 
inter-band transitions between quasi localized states (i.e. 
momentum is not conserved) 46 . Indeed, the method as¬ 
sumes that the dipole-matrix element (proportional to 
the oscillator strength) is not k-dependent. Breaking of 
momentum (k-vector) conservation is a fundamental ten¬ 
ant of quantum confinement in NSs. Furthermore, our 
NSs are amorphous and thus lack long range order. In 
addition, the parabolic band edge picture can still be ap¬ 
plied in a NS as discussed in the next section. A devia¬ 
tion from linearity happens at low ahv due to Urbach tail 
regions. For this reason, the estimation of Eg for amor¬ 
phous materials is usually performed for values of a larger 
than 2xl0~ 9 cm 2 (for our samples) where a clear lin¬ 
ear trend in the Tauc plots appears 47 . Furthermore, the 
Eq was calculated from scanning tunnelling spectroscopy 
measurements on Ge QDs to be ~1.4 eV 48 . This value 
agrees with our Eq determined via the Tauc method for a 
similar diameter range. Therefore, we find that the Tauc 
method works well for amorphous NSs and has been ex¬ 
perimentally demonstrated in the literature 27,47 . 


III. THEORY 

The analysis of the experimental data is carried out 
using a variety of theoretical tools. The approach taken 
in this work considers an interface-related correction to 
the EM A. The validity of the EM A has been widely de¬ 
bated in the literature 5 . Generally, it is assumed that 
the parabolic band approximation breaks down at low- 
dimension and thus the EMA cannot be applied. How¬ 
ever, recently it has been shown that a notion of a ‘fuzzy- 
band-structure’ holds for NSs down to ~ 2 nm 49 , and 
that the k • p method does produce the correct symme¬ 
try of a QD 50 . The EMA is a first order approximation 
of the k ■ p generalization that considers transitions only 
at the conduction band minimum (CBM) and valence 
band maximum (VBM). In a NS, transitions mainly oc¬ 
cur near the Brillouin zone center (the T-point), due to 
breaking of the momentum selection rules 6 . This type 
of transition is referred to as pseudodirect 51,52 . From 
a theoretical perspective, pseudodirect behavior allows 
one to ignore phonon events in the formalism. When 
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a 30-band k • p model is considered, spurious solutions 
are found, which can be removed computationally 53 . A 
two-conduction-band k ■ p approach has also been used 54 
to describe conduction band hybridization, and a two- 
band approximation 6 demonstrates good agreement with 
experiment 25 . A justification of the EMA is given in de¬ 
tail in Refs. 5 and 6 . DFT models including interface 
corrections have been employed in the literature 21,36 , and 
produce a similar dimensional dependence as the method 
presented here 35 . However, we opt for the present ap¬ 
proach because of the clarity in the physics at the inter¬ 
face. Also, the DFT approaches use a different interface 
passivation from the actual experiment one. 

First, we employ a new theoretical formalism using a 
spatially dependent effective mass (SPDEM), as detailed 
in Ref. 30 and formulated in Refs. 55 and 56. The 
important feature here is the introduction of an inverse 
characteristic length scale: y 2 = 2D -2 , where D is the 
QD diameter in nanometres. 7 is not tunable, but was 
derived in Ref. 30 and modifies the dispersion relation 
of the charge carriers as described in Sec. V. 7 acts 
as a coupling parameter between the momentum of the 
charge carriers and the confinement potential. In this 
formalism, the confinement potential, V G , was given by: 


V C , 


,1 


—V 0 exp 



ln 2 (l + 7 i?y) 


(1) 


where r* is the particle position in the i = x,y,z direc¬ 
tion, and V 0 is defined as the energy difference between 
the QD and the matrix material at the CBM or the VBM 
for an electron (V 0 ^ e ) or hole (V 0 ^), respectively. Using a 
point canonical transformation (PCT), we demonstrated 
that the SPDEM is directly related to the confinement 
potential. This effect is analogous to the relationship be¬ 
tween the crystal potential and the bare electron mass 30 . 
Therefore, the SPDEM formalism is well adapted to de¬ 
scribe the effect of V G on the effective mass and yields 
the dispersion relation: 


E g (D) = E g { 00 ) + -j=-^ 




( 2 ) 


where m* e , and m* h are the bulk electron and hole ef¬ 
fective mass, respectively. 

The second model used in our analysis is the effec¬ 
tive mass approximation (EMA), as detailed in Ref. 6 . 
This model describes electron-hole confinement condi¬ 
tions within an infinite confinement potential context, 
according to the relation: 


7.8 


E g (D) = E g ( 00 ) + eV • nm 2 . 


(3) 


We note that EMA presented here is used to describe 
only the effect of reduced dimension on the electron 
and hole and includes no additional effects. However, 
we discovered a physically significant phenomenologi¬ 
cal correction 43 to both the SPDEM and EMA models, 


which was applied to obtain agreement between the two 
theories, given by: 

K D ) = HoCspdemD (l + aD 2+ bD+c ^ ; (4) 

/i(D) = HoCEMA + aD^+bD+c) > (^) 

where a=0.047 nm -2 , b=0.160 nm -1 , and c=-0.035 are 
identical for both renormalized masses, fl{D) and /x(D), 
and were determined in Ref. 43. The parameters ‘a’, ‘b’, 
and l c’ were determined from PL measurements on MBE 
Ge QDs 43 . Ref. 57 noted that these QDs experience 
strong exciton localization and thus there is no Stokes 
shift between the absorption and emission spectrum. Ac¬ 
cording to the theory of QC, localization in position space 
breaks the momentum selection rules and thus phonon 
coupling is not required for indirect gap materials 5 . In 
Refs. 5 and 6 we noted that a measured Stokes shift 
was a consequence of interfacial phononic modes (e.g. 
vibrons). Strong breaking of the momentum selection 
rules (or time-reversibility) denoted as ‘strong confine¬ 
ment’ is included in both the EMA 6 and SPDEM 30 mod¬ 
els. We showed previously 6 that amorphous NSs experi¬ 
ence strong confinement in the electron and hole states. 

In our model, the interface potential influences the ef¬ 
fective mass, as noted above and discussed in Sec. V. 
The parameters cspdem in [nm -1 ] and cema renor¬ 
malize the bulk reduced mass, /z 0 , and will be discussed 
in Sec. IV. Therefore, we conclude that the parameters 
‘a’, ‘b’, and ‘c’ are valid for both absorption and emission 
spectrum, while the influence of the interface is contained 
in the parameter ‘e’. 

The parameters ‘a’, ‘b’, and ‘c’ are the same for all Ge 
QD samples and their physical significance is discussed 
further in Sec. V. Our aim in this work is to understand 
the influence of the interface on the EM. This influence is 
contained solely in the parameter ‘e’, which is introduced 
as correction to the EMA and SPDEM models. This pa¬ 
rameter gives us access to the macroscopic features of 
the interface that would otherwise be inaccessible within 
the standard EMA. Therefore, when we use Eqs. (4) 
and (5) with Eqs. (2) and (3), respectively, we obtain 
perfect agreement in the dispersion relations 43 allowing 
us to isolate essential differences in the physics between 
the two models, as described Sec. V. It is important to 
note that we are presenting a phenomenological analysis 
and not presenting these models as an alternative to pre¬ 
viously existing theoretical models. Our intention is to 
compare the results of two models across different fabri¬ 
cation methods to extract physically relevant information 
correlated with the variation in our parameter set. 


A. Methodology 

For the sputter and PECVD Ge QDs, we use Eqs. (4) 
and (5) in the SPDEM-/7(D) and EMA-/j(P) (see Ref. 
43) models to determine the relationship between the in¬ 
terface potential and EM. For all QD samples the term 













4 



-■l-.-r;.-. -= = = 1 £ - ,^- 1 '*---' _I = - i^JGfc - 

36 34 32 30 28 

Binding Energy [eV] 


FIG. 1. XPS measurement of the Ge 3d peak from the as- 
deposited Ge-rich layer (SiGeO) formed by sputtering. The 
relative intensity of the Ge-Ge peak with respect to its oxide 
states is shown by the respective fits. Color online. 


(l + (aD 2 + bD + c) in Eqs. (4) and (5) is the same, 
thus ensuring that all samples have the same dimensional 
dependence as found in Ref. 43. Therefore, we are left 
to determine only V 0>e , V 0 ,h, and the parameter e. 

First, we consider V 0 , e and V 0t h- From x-ray photoe¬ 
mission spectroscopy (XPS) in sputtered Ge QDs, a Ge- 
Ge signal was measured in the as-deposited film, see Fig. 
1. This signal indicated that Ge QDs formed upon initial 
film deposition. From the XPS results, no Ge-Si nor Si-Si 
signal are observed indicating that all of the Ge is either 
within the QD or in an oxide state, and all of the Si is 
in the Si 02 state 9 . Therefore, sputter thin film depo¬ 
sition, under our given experimental conditions, can be 
modelled approximately by the random mixture model 
(RMM) 44,58 . This model supposes that the film is de¬ 
scribed by the formation of a mixture of Ge, GeO^., and 
Ge 02 centres 59 . A representation of the RMM is given 
in Fig. 2. In the top part of Fig. 2, we indicate a region 
of Ge surrounded by layer of GeO^ and a thin layer of 
Ge 02 - Upon annealing, the Ge QD diameter and the sto¬ 
ichiometric oxide concentration increase. This is a result 
of Ge diffusing from the Ge-rich GeO^ region into the 
Ge QD region 60 , which is indicated in the bottom part 
of Fig. 2. Simultaneously, it is possible that O diffuses 
toward the Ge 02 region. This observation implies that 
the interface between the Ge QD and the matrix mate¬ 
rial is comprised of a large concentration of GeO^ states. 
As a first approximation, we consider a Ge-GeO interface 
with a potential of approximately Vo,e = 0.6 eV and Vo t h 
= 1.8 eV 61 . 

In the case of PECVD QDs, a different behaviour was 
observed. Here the QDs are more accurately modelled 
by the random bonding model (RBM), under our given 
experimental conditions 58 , see Fig. 3. In the RBM all 
of the Ge is assumed to be in one of its oxide states 63 . 
Raman measurements (see Ref. 8 ) indicate that upon 



FIG. 2. Schematic representation of the ideal random mix¬ 
ture model. The top figure depicts the as deposited Ge-rich 
layer where Ge QDs are seen before annealing. Ge QDs are 
surrounded by a sub-oxide shell and pure Ge02- After an¬ 
nealing, the Ge QD diameter increases as stoichiometric Ge02 
forms and the GeO x layer thickness decreases. The diffusion 
process is indicated by the arrows in the middle part of the 
figure. Molecular structures were generated using VESTA 62 
and represent only a cross section of the QD structure. Color 
online. 

deposition the film is primarily in a Ge sub-oxide state, 
i.e. there is no clear observation of a Ge-Ge signal, as 
represented in the top part of Fig. 3. After annealing, 
GeOa, will phase separate into Ge QDs and Ge 02 medi¬ 
ated by Ge diffusion, as represented in the bottom part 
of Fig. 3. As in the case of sputter QDs, it is possi¬ 
ble that O undergoes diffusion as well. This behaviour 
implies that the interface is mainly composed of Ge 02 - 
Though there must still exist a small concentration (as 
compared to sputtered Ge QDs) of GeO^ right at the in¬ 
terface, which will lower the interface potential slightly 
from the ideal Ge 02 value. Therefore, we assume that 
the interface potential is close to the value of Ge-Ge 02 , 
which is Vo, e = 1-2 eV and Vo,h = 3.6 eV 64 . 

The final step is to determine e for both the SPDEM 
and EM A models. We consider a Tauc plot of the ex¬ 
perimental absorption cross section, <r, for PECVD and 
sputtered Ge QDs from Refs. 8 and 9, respectively, a is 
modelled via the relation: 

c = {hu - E G {D)f ; (6) 

where B* is proportional to the the oscillator strength 
and was determined experimentally, and hv is the pho¬ 
ton energy. Eq(D) is given by either SPDEM- fl(D) or 
EMA-^(D), in Eq. ( 6 ). Eq{D) and D was determined 
experimentally in Refs. 8 and 9, see Sec. II. In the case of 
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FIG. 3. Schematic representation of the ideal random bonding 
model. The top figure depicts the as deposited Ge-rich layer 
comprised of sub-stoichiometric Ge02, GeCU, where x <1.5. 
After annealing, phase separation occurs and Ge QDs form 
surrounded by Ge02, with a small concentration of GeCU 
right at the interface. The diffusion process is indicated by the 
arrows in the middle part of the figure. Molecular structures 
were generated using VESTA 62 and represent only a cross 
section of the QD structure. Color online. 


EMA-^(D), we fit Eq. (6) via Cema to the experimental 
data for each QD D and Eg{D). For SPDEM-/I(.D), as 
we noted above, the interface potential was approximated 
from our discussion of the RMM, or RBM for the sput¬ 
ter, or PECVD samples, respectively. Since these models 
are only a first approximation we allowed V 0 , e and V a< h 
to vary within 25% of their proposed values, and then we 
fitted Eq. (6) to the data by varying espDEM and V oe (h\ 
for each experimental QD D and Eg{D). 


IV. RESULTS 

The results obtained from fitting SPDEM-/2(D) to the 
sputtered and PECVD samples are shown in Figs. 4 and 
5, respectively. The deviation from linearity in the ex¬ 
perimental data at low ahv is due to Urbach tail regions, 
i.e. defect states near the band edge, which are common 
to amorphous structures as in this work 9 . In PECVD 
samples (Fig. 5) there is a deviation from linearity at 
high <jhv due to a partial crystalline transition in less 
than 30% of the QDs after annealing at 800 °C 8 . When 
the fit was constrained to the linear region, a less than 
1 % difference was found in our final fitting parameters 
compared to fitting over the entire data set. The fitted 
value for cspdem is listed for each QD D in Figs. 4 and 



FIG. 4. Tauc plot of the absorption cross section with corre¬ 
sponding linear fit for sputtered Ge QDs from Ref. 9. The 
diameter for each sample is given in the legend, with the fit¬ 
ting parameter e from Eq. (4). The fitted interface energy 
from Eq. (2) is labelled in the plot. 



FIG. 5. Tauc plot of the absorption cross section with cor¬ 
responding linear fit for PECVD Ge QDs from Ref. 8. The 
diameter for each sample is given in the legend, with the fit¬ 
ting parameter e from Eq. (4). The fitted interface energy 
from Eq. (2) is labelled in the plot. 


5, and the resultant interface potential (also shown), was 
found to be notably identical (within 1%) for each QD 
D. espDEM for each QD D are the same within experi¬ 
mental error, therefore, in our formulation of Eq{D) we 
use the average value, (e)=0.181, and 0.090 for sputter, 
and PECVD, respectively. Likewise, Cema varies only 
within experimental error. The results of our analysis 
are summarized in Table I. Table I also lists the results 
from the EMA fit and the values for MBE Ge QDs. 

In Fig. 6, we plot the experimental data for PECVD 8 
and sputter 9 Ge QDs along with the SPDEM-/I(D) curve. 
The interface potential, Cema , and Eq{ oo) are listed in 
Fig. 6 for reference. In Sec. Ill, we mentioned that the 
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TABLE I. Renormalized mass parameter, e, for SPDEM (Eq. 
(4)) and EMA (Eq. (5)), and interface potential from Eq. (2) 
for sputter, PECVD, and MBE Ge QDs. 




Sputter 


2 nm 

2.5 nm 3 nm 4 nm 

<e> 

e-EMA 

0.528 

0.597 0.583 0.499 

0.552 

esPDEM (nm -1 ) 

0.173 

0.195 0.191 0.163 

0.181 


v J 

V °i 

'“=0.92 eV V/%= 2.80 eV 



PECVD 



3.5 nm 

4.4 nm 8.4 nm 

(e) 

e-EMA 

0.203 

0.268 0.288 

0.253 

espDEM (nm -1 ) 

0.072 

0.095 0.102 

0.090 


vl 

'“=1.09 eV V/“=3.28 eV 



MBE 


Gem a 


0.378 


espDEM (nm -1 ) 


0.059 



V a 

lje =0.20 eV V o ,h=0.66 eV 




FIG. 6. Variation in the gap energy as a function of QD 
diameter for PECVD 8 and sputter 9 Ge QDs. The theoretical 
curves are from Eq. (2) with the correction factor Eq. (4). 
The MBE curve is taken from Ref. 43 for reference. The 
bulk gap energy Eg{ oo), interface energy, and e is given for 
each sample from Table I. The MBE curve uses the Ge/a-Si 
interface energy 43 . 


SPDEM-/I(D) and EMA-/i(D) models agree with each 
other. The parameters a, b , and c in Eqs. (4) and (5) 
were found by analysing MBE grown Ge QDs 43 . There is 
good agreement between the experimental data and our 
calculation, noting that the error bars in the D of the 
PECVD samples are on the order of 2 nm 8 . 

The Eg for the PECVD grown samples is higher than 
the sputter samples. The reason for the difference is un¬ 
derstood within our models. From the SPDEM-/i(D) 
model, we find a larger interface potential in PECVD 
samples compared to the sputter samples. The larger 
interface potential is correlated with a larger reduction 


in 6emAi which increases the confinement energy. MBE 
Ge QDs were fabricated with an a-Si matrix, which has 
a lower interface potential than both PECVD and sput¬ 
ter Ge QDs, see Fig. 6 . However, the Eq is compa¬ 
rable to the other two samples with a value of cema 
in between the PECVD and sputter values. It is diffi¬ 
cult to make a full comparison between the sputter and 
PECVD samples versus the MBE samples, because the 
MBE samples were fabricated under very different con¬ 
ditions. However, we can consider a few essential differ¬ 
ences. In the case of the MBE QDs there are no O-defect 
states 43,65 , because of the unique fabrication method. O- 
defect states at the interface increase the spread of the 
wave-function, thus lowering the confinement energy 12 . 
Furthermore, O interface states can pin the charge carri¬ 
ers, again lowering the confinement energy 5,6 . The com¬ 
bination of these phenomena lowers the confinement en¬ 
ergy in PECVD and sputter samples, in agreement with 
theoretical predictions 12 , and it becomes comparable to 
that in MBE samples (see Fig. 6 ). Therefore, we under¬ 
stand that the O interface states are the source of the 
difference between the sputter and PECVD samples ver¬ 
sus the MBE samples. Nonetheless, we see that for a 
higher interface energy a larger reduction in the EM is 
obtained. 


V. DISCUSSION 

A closer examination of the results shown in Table I 
in the context of our models allows us to understand the 
nature of the confinement mechanism. As mentioned in 
Sec. Ill, the SPDEM theory models a modification in the 
EM due to the interface potential. In Ref. 30, we derive 
the envelope function within the SPDEM formalism: 

^i(D) = ex P 9^2 + 7 i r i)j ; (7) 

where a ~ (V ),) -1 / 2 is the Gaussian strength. In Fig. 
7, we plot the probability amplitude of Eq. (7) for a 
6 nm QD as a function of the interface potential and 
particle position. As the interface potential increases, 
the probability amplitude is increasingly localized, which 
implies increasing de-localization in momentum space. A 
Fourier transformation of Eq. (7) demonstrated that the 
SPDEM envelope function exhibits an increased spread 
in momentum space compared to a Gaussian envelope 
function 30 . The interface potential alters the SPDEM 
envelope function because of the introduction of 7 , which 
modified the EM through the relation 30 : 


Therefore, from ( espDEM ) _1 we can find the length 
scale at which the EM begins to change, which is 16.95, 
11.11, and 5.53 nm in MBE, PECVD, and sputter sam¬ 
ples, respectively. At the stated dimensions the SPDEM 
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Probability Amplitude: 6 nm Quantum Dot 
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FIG. 7. Probability amplitude of Eq. (7) taken from Ref. 30 
for a 6 nm Ge QD as a function of particle position and in¬ 
terface potential. The figure depicts increasing localization of 
the probability amplitude as a function of increasing interface 
potential. Color online. 


energy begins to deviate from the experimental data and 
&spdem corrects for this behaviour. In Fig. 8 we plot 
p,(D)/n 0 and n(D)/fj, 0 from Eqs. (4) and (5), respec¬ 
tively, where the role of s-spdem can be seen. Below di¬ 
ameters of around 16 and 11 nm the MBE and PECVD 
fl(D)/pi 0 curves are less than 1 indicating that jl(D) in¬ 
creases Eg{D) below these dimensions. The difference 
between the cross-over dimension depicted in Fig. 8 and 
the value obtained from ( esPDEM)~ 4 is because of exper¬ 
imental error. We note that in Fig. 8 the curve 

for the sputter sample does not fall below 1. In Ref. 43, 
we determined that above 6.2 nm the dispersion relation 
is dominated by the SPDEM-/i(i4) relation. Therefore, 
the sputter {cspdem) 1 value of 5.53 nm indicates that 
the SPDEM mechanism is suppressed in these samples. 
The different behaviour found for the sputter samples 
is a consequence of their interface quality (Secs. Ill A 
and IV). The lower interface potential does not create a 
strong confinement condition for the charge carriers, see 
Fig. 7. Additionally, the sub-oxide interface states pro¬ 
mote de-localization of the wave-function in real space 12 . 

The values obtained for ( csdpem)~ 1 also agree with 
the results of Fig. 6. First, we note that the value ob¬ 
tained for the MBE samples of ~ 16 nm is on the order 
of the Bohr radius for Ge. This result is what we expect 
in these samples, because of the high quality interface. 
At 16 nm, the charge carriers in the MBE samples begin 
to experience the influence of the confinement potential 
and thus their energy is modified, as predicted by QC. 
While in the PECVD and sputter samples, the onset of 
QC is suppressed from the expected Bohr radius value. 
Therefore, the PECVD and sputter samples experience 
the effects of QC over a smaller diameter range compared 
to MBE, and thus their Eg shift is not as great. 

The relationship between pure QC effects and a modi¬ 
fied EM are understood by comparing the SPDEM-/I(£)) 



FIG. 8. Plot of Eqs. (4) and (5) normalized to the bulk 
reduced effective mass, fi 0 , for the different Ge QD samples 
as a function of QD diameter. The SPDEM-/i(D) (Eq. (4)) 
curves are indicated by an arrow. The vertical lines indicate 
the dimension where the behaviour changes from SPDEM to 
EMA. 


results with respect to the EMA-/x(D) results. In a pre¬ 
vious work on Ge QDs buried in a-Si 43 we found that 
above «6.2 nm the Eq ~ D -1 , due to a modification 
in the EM, while below rs 4.6 nm the Eq ~ D~ 2 , due 
to pure QC effects. These parameters come from the 
term ^1 + (aD 2 + bD + c) ^ in Eqs. (4) and (5), and 

represent a universal dimensional dependence. The pa¬ 
rameter e represents a correction to these ideal length 
scales due to additional effects not included in the the¬ 
ory, which in our analysis include the O-interface states. 
If we consider the ratio Cema /egp deMi then we obtain 
3.05 and 2.81 nm for sputter and PECVD samples, re¬ 
spectively. This length scale signifies when the confine¬ 
ment behaviour changes from the SPDEM to the EMA 
model and is depicted in Fig. 8 by the vertical line at 2.9 
nm. Within experimental error, we thus find that pure 
QC effects dominate at « 2.9 nm, which is suppressed 
from the ideal case (between 6.2 and 4.6 nm) due to the 
interface states. On the other hand, the same ratio in 
MBE samples gives 6.4 nm (see vertical line in Fig. 8), 
which is in agreement with the ideal value (6.2 nm), be¬ 
cause there are no O interface states in this sample. Fi¬ 
nally, this is in agreement with previous results obtained 
in Ref. 43, where the MBE value of 6 nm was found to 
be the dimension that produces the most efficient lumi¬ 
nescence. From our work here we understand this result 
on the basis that this is the dimension at which the be¬ 
haviour changes from SPDEM to EMA and thus carriers 
are strongly confined. 

We illustrate collectively all of the phenomena de¬ 
scribed in this manuscript in Fig. 9. At the top of 
the figure we recall the dimensional dependence of the 
SPDEM-/i(D) and EMA -p>(D) models. Fig 9 is divided 
into two parts. In the ‘pink’ region of the figure above ~ 6 























nm, the dominant model is SPDEM-/i(-D), which models 
a change in Eq in conjunction with a change in the EM. 
The ‘grey’ region below ~ 4 nm is modelled according 
to the EMA-/i(D) where pure QC effects dominate. The 
‘pink’ and ‘grey’ regions were determined from the results 
of Ref. 43 and come from the universal dimensional pa¬ 
rameters a , and b in the term ^1 + (aD 2 + bD + c) ^ . 

From the parameter (cspdem) -1 , we determined the 
onset of the SPDEM behaviour at ~ 16, 11, and 5 nm 
in the MBE, PECVD, and sputter samples, respectively. 
The change in energy due to SPDEM is represented by 
the solid curved line in Fig. 9 for the three samples (not 
drawn to scale). The MBE sample’s value is in agree¬ 
ment with the Bohr radius for Ge and indicates when 
the charge carriers first ‘feel’ the confinement potential. 
In the case of PECVD and sputter samples this value is 
suppressed due to the O states at the interface, which 
cause the wave-function to leak into the matrix and thus 
lowers the confinement energy, as described above. The 
PECVD case has a higher value than the sputter case due 
to the higher confinement potential and thus experience 
a larger change in Eg, due to an increased reduction in 
the EM. 

The transition point from SPDEM to EMA is on the 
order of 5 nm according to the transition from the ‘pink’ 
to the ‘grey’ region in Fig. 9. From the ratio of cema 
to s-spdem we determined this transition for each sam¬ 
ple. The change in energy due to EMA is represented by 
the dashed curved line in Fig. 9 for the three samples 
(not drawn to scale). In the case of MBE sample, the 
transition is in agreement with the calculated value at 6 
nm, because of the lack of O at the interface. While in 
PECVD and sputter samples they both change behaviour 
at around 3 nm. The suppression here again is due to O 
interface states, as described above. 


VI. CONCLUSION 

Using a previously developed theoretical formalism as 
an analytical tool we have been able to understand the 
relationship between the EM and the interface potential. 
We found that when the interface potential was large 
there exists a strong reduction in the EM. This condition 
was clearly demonstrated when comparing PECVD and 
sputter samples, where the PECVD case has a larger 
interface potential. However, if the sample contains O 
states at the interface, then the confinement energy was 
reduced due to these states. This result was noted by 
comparing MBE samples with no O interface states to 
the PECVD and sputter samples. We determined the 
length scale where the samples begin to experience a 
modification in Eg due to a change in the EM, which 
was ~ 16, 11, and 5 nm for MBE, PECVD, and sput¬ 
ter samples, respectively. The lower values in PECVD 
and sputter samples compared to MBE samples were un¬ 
derstood in terms of the O interface states. We found 


the length scale where the QC mechanism changes from 
SPDEM to EMA, which was 6, and 3 nm in MBE and 
PECVD/sputter samples, respectively. All of these re¬ 
sults are summarised graphically in Fig. 9. These re¬ 
sults indicate an additional carrier confinement mecha¬ 
nism in NSs, which can be exploited for device fabrica¬ 
tion through interface engineering. From this work, we 
have found a dynamic relationship between the magni¬ 
tude of the interface potential and the EM with respect 
to the chemical composition of the interface. 
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